Competition of synchronization domains in arrays of chaotic homoclinic systems. 
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We investigate the response of an open chain of bidirectionally coupled chaotic homoclinic systems 
to external periodic stimuli. When one end of the chain is driven by a periodic signal, the system 
propagates a phase synchronization state in a certain range of coupling strengths and external 
frequencies. When two simultaneous forcings are applied at different points of the array, a rich 
phenomenology of stable competitive states is observed, including temporal alternation and spatial 
coexistence synchronization domains. 
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Synchronization of chaos refers to a process wherein 
two (or many) chaotic systems adjust a given property 
of their motion to a common behavior ll. The emer- 
gence of synchronized features has been investigated in 
nature 0, in controlled laboratory experiments 0, Q, 
and the interest has moved toward the characterization 
of synchronization phenomena in spatially extended sys- 
tems, such as large populations of coupled chaotic units 
and neural networks [5|, globally or locally map lattices 
and continuous space extended systems 0). 

In this Letter, we consider a one dimensional chain of 
sites, each one undergoing a local homoclinic chaotic dy- 
namics, interacting via a bidirectional nearest neighbor 
coupling. Homoclinic chaos consists of a train of nearly 
identical spikes separated by erratic inter-spikes inter- 
vals (ISI). In phase space, this motion corresponds to the 
passage through a saddle focus, where stable manifolds 
collapse and an unstable manifold emerges, with the ex- 
pansion rate larger than the contraction one The 
saddle region displays a large susceptibility to an exter- 
nal stimulus, therefore such a chaotic system gets easily 
synchronized to a weak forcing signal. 

The ability of such systems to synchronize to an ex- 
ternal forcing was demonstrated in previous works 
finding that it may constitute a reliable communication 
channel Q robust against noise [lj}. Furthermore, ho- 
moclinic chaos can be self-synchronized by feeding back 
a finite train of its own spikes, via either a delayed feed- 
back ^lj or a low frequency filter 0] j m this latter case 
obtaining a bursting behavior reminiscent of the dynamic 
of neurons in Central Pattern Generators fl3j |. 

However, when passing from a single system to an ar- 
ray, a relevant problem emerges related to the ability of 
the array to respond to external periodic perturbations 
localized at one end site, yielding synchronized patterns. 
The issue we are addressing is relevant for biological or 
artificial communication networks. 

For convenience, we refer to a chain of dynamical units, 
each one represented by a 6- variable system modeling ho- 
moclinic chaos in a single mode CO2 laser with feedback 



[lij . The extension of the model to an array is: 

x\ = k x\(x 2 - 1 — ki sin 2 x l 6 ), 
±2 = -71^2 _ 2k x\x 2 + gx\ + x\+p, 
±1 = -jix l 3 + gx\ + 4 + p, 
x\ — -72X4 + zx\ + gx\ + zp, 
x\ = -722:5 + zx\ + gx\ + zp, 

x" 



(1) 



= -f3[x l 6 -b + r(- 



+ e{x\~ 1 +x\ +1 - 2< x\ >))] 



Here the index i denotes the i th site position (i = 1, .., N), 
and dots denote temporal derivatives. For each site, xi 
represents the laser intensity, xi the population inversion 
between the two levels resonant with the radiation field, 
and Xe the feedback voltage which controls the cavity 
losses. The auxiliary variables 2:3,2:4 and 2:5 account for 
molecular exchanges between the two resonant levels and 
the other rotational levels of the same vibrational band. 
We consider identical units; as for the parameters, their 
physical meaning has been already discussed (l^j . Their 
values are: k Q = 28.5714, ki = 4.5556, 71 = 10.0643, 
72 = 1.0643, g = 0.05, p = 0.016, z = 10, = 0.4286, 
a = 32.8767, r = 160, b = 0.1032. 

The coupling on each site is realized by adding to the 
xq equation a function of the intensity (2:1) of the neigh- 
boring oscillators. The term < x\ > represents the av- 
erage value of the x\ variable, calculated as a moving 
average over the whole evolution time. The coupling 
strength e > is our control parameter. The system is 
integrated by means of a standard fourth-order Runge- 
Kutta method with open boundary conditions. 

We first study the emergence of synchronization in the 
absence of external stimuli, as the coupling strength e 
increases. Due to the coupling, a spike on one site induces 
the onset of a spike in the neighbor sites as discussed in 
Ref. [3. 

In Fig. ^a)-(b) we show the transition from unsyn- 
chronized to synchronized regimes by a space time repre- 
sentation of the array. A detection threshold isolates the 
spikes getting rid of the chaotic small inter-spike back- 
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FIG. 1: Space-time representation of spikes positions for e=0.05 
(a) and e=0.2 (b). (c) Average defect density vs. e, for different 
chain lengths: N=10 (□) , 40 (o), 80(*). 



ground, thus we plot only the spike positions as black 
dots. The transition to phase synchronization is antici- 
pated by regimes where clusters of oscillators spike quasi- 
simultaneously [l^. Clusters are delimited by "phase 
slips" or defects, easily seen as holes in the space time 
fabric. More precisely, we introduce a phase measure 
4> l (t) for a time interval t between two successive spikes 
of the same site, occurring at T % k _ x , t^, by linear interpo- 
lation 0: 

= 2*4=^-. (2) 

T k T fc-1 

A defect appears as a 2ir "phase slip" in the difference 
between the phases of two adjacent sites. Notice that this 
mutual referencing is the natural extension of measuring 
the regularity of a sequence against an external clock, 
whenever there is no external clock, but the time evolu- 
tion of a site compares with the nearest neighbor sites. 
Hence "phase synchronization" denotes a connected line 
from left to right, not broken by defects. This defini- 
tion of phase synchronization does not imply equal time 
occurrence, thus the unbroken lines are not isochronous 
as can be seen in Fig. ^ (b). The cluster size increases 
with e extending eventually to the whole system (Fig. 
IHb)). The route to phase synchronization can be char- 
acterized by the defect density, that is, the number of 
defects per site. In Fig. QJc) we plot the average defect 
density as function of e, calculated for a long evolution 
time (T = 3 x 10 5 ). Full phase synchronization is reached 
when the defect density falls below one defect per site. 
The defect statistics has been studied for several chain 
lengths; we find that above iV=30 there are no apprecia- 
ble size-dependent effects. Once phase synchronization is 
established, a further increasing of e reduces the natural 
frequency of the system ui (e) = < /g 2 / > ( e ) ■ We will show 
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FIG. 2: Response of the N=40 chain with e = 0.13 and 
cj o (e)=0.02 to an external periodic forcing: (a) w=0.015 , (b) 
u>— 0.042. (c) Penetration depth vs. u for different coupling 
strengths e=: 0.12 (*), 0.15 (o), 0.2 (□), 0.25 (o). 



that this slowing down affects the capability of the array 
to synchronize to an external signal. 

For this purpose, we explore the response of the sys- 
tem to an external periodic stimulus applied to the first 
site of the chain. Precisely, we periodically modulate 
the parameter b at the site i = 1 as b\ — b a (t) = 
b a (\ + Asm(ujt)). From previous work, we know that 
this driving can induce a phase synchronization on a sin- 
gle oscillator Q ; here we explore the ability of the system 
to transmit the periodic signal through the chain. 

The modulation amplitude A does not affect the re- 
sults, provided that it is sufficient to synchronize the 
i = l site. Therefore, we do not loose generality by fixing 
a constant value A = 0.3. 

We will consider that the signal has been successfully 
transmitted through the system when after a finite time 
the last element of the chain spikes with the same period 
of the external forcing, without defects. In Fig. |21(a)-(b) 
we give examples of partial signal transmission. If uj is 
too small (Fig. [21 (a)) or too large (Fig. [21 (b)) with re- 
spect to the natural frequency of the system (w o (e)=0.02 
in the figure), then only partial transmission is achieved. 

We explore the (e, uj) range over which transmission 
propagates over the whole chain. In the low frequency 
limit, we find that u) < uj (e) is not able to globally syn- 
chronize the chain. Independently of N, as uj d is larger 
than uj, the last sites tend to spike spontaneously between 
two consecutive periods of the external driver before the 
synchronization propagates to them, and therefore syn- 
chronization is lost (Fig. [21(a)). 

When uj > 2ui (e), the first N p sites synchronize with 
the driving frequency, but beyond N p a line of defects 
restores the natural oscillation regime (Fig. [21(b)). This 
"penetration depth" N p for synchronization is invariant 
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FIG. 3: Curves delimiting the (e, ui) range for synchronized 
transmission in the arrays with lengths: N=10 (dotted lines), 
N=20 (dotted-dashed lines), N=40 (dashed lines) and N=80 
(solid lines). The black stars indicate the average spiking 
frequency uj (e) of a site of the array in the absence of external 
perturbation. Inset: zoom of the area where the two external 
frequencies are selected for studying the spatial competition 
between synchronization domains. 
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FIG. 4: Competition between spatial synchronization regimes 
induced by external forcing: (a) ivi = 0.02, u)2 = 0.021, e = 
0.13, (b) u/j = 0.038, uj 2 = 0.042, e = 0.12, (c) wi = 0.04, w 2 = 
0.0405,e = 0.11 



as we change the whole array length. In Fig. 121(c) we plot 
the penetration depth vs. the forcing frequency for differ- 
ent values of e. If for a point (e, ui) the penetration depth 
is N p , then one would observe complete synchronization 
only for sizes N < N p , while incomplete synchronization 
would unavoidably take place for N > N p . As a result, 
for a given N, only a limited range of external frequencies 
can be transmitted over the whole chain. 

In Fig. |2| we plot the boundaries of the transmission 
band as a function of e and u>, for several chain lengths. 
The region inside the curves contains all the (e, to) points 
for which global transmission is allowed. It can be seen 
that for each e, the transmission band extends from io 
(black stars) to approximately 2lo q . Notice that the sys- 
tem starts to transmit for coupling strengths above the 
ones leading to intrinsic synchronization (approximately 
e > 0.11). For weaker couplings, the presence of defects 
breaks the continuity, while for e > 0.35 the homoclinic 
dynamics is destroyed. The left boundary of the trans- 
mission range refers to perfect transmission of the u pe- 
riod up to the end of the chain. If one is only interested in 
the transmission of the average frequency, this boundary 
is slightly smeared out. 

Now we have sufficient background to address the main 
question of how two different frequencies applied at the 
far ends of the chain compete in generating two separate 
spatial patterns of synchronization. The temporal com- 
petition between different synchronization states was re- 



cently investigated theoretically [TjJ and experimentally 
[lij in the context of a single chaotic system forced by two 
external frequencies, finding competitive behaviors as al- 
ternations of synchronism to several frequencies {uj\ , u>2 
or a combination of the two) . Here we go further address- 
ing the problem of spatial competition between synchro- 
nization regions, which is preliminary to controlling the 
dynamics of an extended system as well as to studying 
the response of neural assemblies to competing external 
perturbations. 

To answer this question, we apply to the first (i = 1) 
and last (i — N) site two periodic perturbations with 
frequencies lu± and u)%, respectively. For simplicity, we 
will take always ui < uj\ < u>2, so that N p (ui) > N p (u/2) 
(see Fig. 0(c)). 

The emerging competition scenario can be described 
with reference to Fig. El For N p (u)i), Npfa) > N, both 
frequencies synchronize over the whole chain. However, 
after a suitable transient time, the whole system synchro- 
nizes to the larger frequency u>2 , with the only exception 
of the site i = 1 (Fig. 0] (a)). This kind of "winner- 
takes-all" behavior is the consequence of the extended 
character of the system, and is at variance with the sin- 
gle oscillator behavior for both forcing frequencies inside 
the Arnold tongue. In fact, in Ref. [lj| the entrain- 
ment takes place at the frequency closer to the natural 
frequency lo . 

For N p (uji) > N,N p (lu 2 ) < N, only the smaller fre- 
quency synchronizes over the whole chain, while the 
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FIG. 5: (a) Time evolution of co(t) for i=l and i=40 (dashed lines) 
and for the site i=22 (solid line) located on the domain boundary 
for the case of Fig. ||] (b). (b) Normalized locking time at the 
frequencies u)\ (*) and W2 (o) vs. site index for the whole chain. 



larger frequency is limited to the N p (ui2) sites closest to 
i = N. In this situation, we find that permanent syn- 
chronization domains for oj\ and u>2 are established, with 
an irregular domain wall (Fig0] (b)). If we increase N , 
the loi domain is always confined to the last N p (w2) sites, 
independently of the total length of the chain as well as 
of the value of u>\. The domains are stable, as checked 
with a very long integration time (T > 3 x 10 7 ). For the 
case of Fig 01 (b), we plot the instantaneous frequency 
for i = l,i = 40 and i = 22 (Fig. (a)). It can 
be observed that the site i = 22 located on the domain 
boundary locks alternatively to lo\ and loi. The lock- 
ing periods are interrupted by defects. In Fig 0(b), we 
plot the normalized locking time of all sites to u>i and 
u)2, respectively. The transition is smooth spacewise and 
the boundary layer has a width of approximately 6 sites, 
independently of the chain length. 

Finally, for N p (wi_) , N p (u)2) < N, neither of the two 
frequencies stabilizes a synchronized pattern over the 
whole chain. In this case, we observe alternation be- 
tween synchronization patterns with frequencies u>i and 
u>2, with intervals of asynchrony filled with defects, as 
shown in Fig (c) . The duration of the synchrony and 
asynchrony intervals is irregular. This competitive be- 
havior persists in time. As a consequence, the competi- 
tion between the two frequencies has here a cooperative 
effect, insofar as it enhances the ability of each single en- 
trainment process to reach global synchronization over 
finite time slots. 

In summary, we have studied the response of a chain of 
nearest neighbor coupled homoclinic oscillators to peri- 
odic stimuli. The array can propagate a synchronization 



state in a range of couplings and external frequencies 
(e,aj). When two simultaneous forcings are applied at 
different points of the array, a rich phenomenology of 
stable competitive states is observed. The features and 
stability of these states depend on the intrinsic dynamics 
of the system independently of the chain size. 

The Authors are indebted to R. Meucci for fruitful dis- 
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CT-2000-00158. 
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